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Abstract 



The effect of direct processes on the statistical properties of deterministic 
scattering processes in a chaotic waveguide is examined. The single channel 
Poisson kernel describes well the distribution of S-matrix eigenphases when 
evaluated over an energy interval. When direct processes are transformed 
away, the scattering processes exhibit universal random matrix behavior. The 
effect of chaos on scattering wavefunctions, eigenphases, and time delays is 
discussed. 
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I. INTRODUCTION 



In 1957, Wigner proposed the use of statistical measures to analyze complex nuclear 
scattering data [1]. It was soon found that, in some nuclear scattering data, the spacing 
between scattering resonances was distributed in a manner similar to the spacing of energy 
levels of a Gaussian random Hamiltonian [2]. In the 1960's, extensive work was done to 
develop a systematic theory of the statistical properties of random hermitian matrices [3,4] 
and random unitary matrices [5,6]. The general criterion for constructing these random 
matrices is that they minimize information. In 1979, the appearance of random matrix-like 
behavior in quantum systems was linked to underlying chaos in the classical deterministic 
dynamics of these systems [7-10]. Since then a large body of work has developed linking 
the statistical properties of bounded and open quantum systems to underlying deterministic 
chaos [11,12]. 

Historically there have been two different approaches to describing the statistics of quan- 
tum scattering processes in chaotic systems. One approach [13] begins directly with the 
Hamiltonian and uses it to build the scattering matrix. The other approach [14,15], begins 
directly with the scattering matrix. In both cases, random matrices are used to describe 
scattering processes and the statistical properties of the eigenphases of the scattering matrix 
and partial delay times can be obtained and compared to experiment. There has been con- 
siderable success in recent years to link the results and predictions of these two approaches 
[16]. 

The connection between the statistical properties of scattering processes and underlying 
chaos is not straightforward because scattering events may involve either reactive processes 
or direct processes. Reactive scattering processes are those for which an incident particle 
becomes engaged with the dynamics in the reaction region, and may be delayed there for 
a considerable time. Direct processes are those for which the particle passes through the 
reaction region without becoming significantly engaged in the reaction region. One predic- 
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tion of the random matrix theory (RMT) of scattering is that, if the scattering process is 
truly chaotic, the average value of the S-matrix will be zero. However, if direct processes 
are present this will not happen. 

In this paper, we will study the deterministic scattering of a quantum particle in a two 
dimensional ballistic waveguide which has a classically chaotic cavity formed by a ripple bil- 
liard connected to a single lead at one end (see Figure 1). The ripple billiard is particularly 
well suited to the use of the reaction matrix theory approach to scattering [17], because a 
simple coordinate transformation [18,19] allows us to construct a Hamiltonian matrix and 
thus an eigenvalue equation for the basis states inside the cavity. We will compare the results 
of deterministic scattering from the chaotic cavity to some recent predictions of random ma- 
trix theory as regards the scattering process. The open geometry of the waveguide in Fig. 
(1) ensures that direct processes will play an important role in scattering. We show that the 
contribution of the direct processes to the scattering can be transformed out and the sta- 
tistical properties of the reactive part of the scattering process can be compared to random 
matrix predictions. Our discussion focuses on quantum particles and we will use parameters 
appropriate for electrons in waveguides made with GaAs, for which a number of experi- 
ments have been done [20-22]. However, our results also apply to electro-magnetic waves in 
flat microwave cavities, because the eigenmodes in these cavities satisfy a Schrodinger-like 
equation [23-25]. 

We begin in Section (2), by reviewing the reaction matrix theory of deterministic scat- 
tering in the ballistic waveguide and we study some properties of the cavity basis states. 
In Section (3), we study various statistical properties of the eigenphases of the waveguide 
scattering matrix (S-matrix). We show that, when the cavity dynamics is chaotic, the deter- 
ministic S-matrix eigenphases exhibit level repulsion and their distribution is well described 
by a Poisson kernel. In Section (4), we compare the distribution of partial delay times for 
the deterministic scattering process to the predictions of random matrix theory. Finally, in 
Section (5) we give concluding remarks. 
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II. REACTION MATRIX DESCRIPTION OF SCATTERING 



We will explore the behavior of a particle of mass, m, in a ballistic waveguide as it 
scatters from the ripple cavity shown in Fig. (1). A particle with energy, E, enters the 
cavity from the left along a straight lead which has infinitely hard walls. The particle wave 
is reflected back to the left by an infinitely hard wall located at x — L. The dynamics inside 
the cavity, < x < L, can range from mixed to fully chaotic as the ripple amplitude is varied. 
The Schrodinger equation, which describes propagation of a particle wave, ^(x, y,t), in the 
waveguide at time, t, is given by 



h 2 / d 2 d 2 



*(x,y,t), (2.1) 



2m V dx 2 dy 2 , 

where H is the total Hamiltonian of the particle in the waveguide and % is Planck's constant. 
The potential, V(x, y), has the following properties: V(x, y) = oo for (L<x < oo); V(x, 0) = 
oo for (— oo<x<L); V(x,y = g(x)) = oo for (0 < x < L); and V(x, y = d) = oo for 
(— oo < x < 0); where g{x) = d + a sin(57ra;/L) gives the contour of the ripple, d is 
the average width of the cavity, L is the length, and a is the ripple amplitude. In all 
subsequent sections, we take the particle mass to be the effective mass of an electron in 
GaAs, m = 0.067m e , where m e is the free electron mass. An energy eigenstate, \E), with 
energy, E, satisfies the equation, H\E) = E\E), and evolves as ty(x,y,t) = (x,y\E) e - lEt / h . 
As shown in reference [19] the waveguide energy eigenstates can be expressed in the form 

oo oo 

(x, y \E) = Y,iM^v) + E r « *n(x,y)- (2.2) 

j = l 71=1 

The states, & n (x,y), are the basis states in the lead (x<0), 

$ n (x,y)={x,y\$ n ) = ]J^Xn(x)sm(^^j (2.3) 

(n = l,2,...,oo). These will consist of both propagating and evanescent modes, as we 
will discuss below. The states, (f>j(x, y)=(x, y\4>j) (j = 1,2,..., oo), are the basis states in the 
cavity (0<x<L). In practice, we truncate the number of cavity basis states to some large but 
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finite number, M. The value of M is determined by the accuracy desired for the calculation. 
The coefficients, T n and 7,-, in Eq. (2.2) are defined T n = (& n \E) and lj — (4>j\E)- 

A complete set of basis states, <f>j (x, y) , inside the ripple cavity has been obtained 
by solving the Schrodinger equation in the cavity, using Neumann boundary conditions 
((t 1 )^ = 0) at the entrance (x = 0), and Dirichlet boundary conditions, (<f)j(x,y) = 0) at 
the walls. We obtained the eigenfunctions, <f>j(x,y) and the associated eigenvalues, Xj, using 
the transformation technique discussed in Ref. [19], Sect. V. 

In Fig. (2. a), we show a Poincare surface of section for a classical particle in a closed 
ripple cavity with the same shape as in Fig. 1, and with hard walls. We choose d = 101 A, 
L = 51 lA, and a = I.OA. We have plotted Birkhoff coordinates, p x /p = cos(o;) versus 
x, each time the particle hits the lower wall at point x (a is the angle between the wall 
and the momentum). For these billiard dimensions, the classical phase space contains a 
mixture of regular orbits, nonlinear resonances, and chaotic motion. If we increase the 
ripple amplitude, a, there is a range of values of a for which the classical motion appears to 
become totally chaotic. For the scattering system (with one end of the cavity open) periodic 
orbits and nonlinear resonance regions can't be reached classically by a particle that enters 
from the left, but quantum mechanically tunneling into these regions is possible as we will 
show Section IV (see also [26]). In Figs. (2.b-e), we show Husimi plots of quantum Poincare 
surfaces of section [30] of cavity basis states, <pj(x,y), with eigenvalues 7 983 = 257.1934E 1 !, 
7985 = 257.9339^1, 7989 = 258.6655^1, and 7 990 = 258.9072.E7i, where E l = gg. The 
Husumi plots in Fig. (2.c) and Fig. (2.d) show that these bases states reside on nonlinear 
resonances of the underlying classical phase space. In Section IV, we will show that these 
bases states give the primary contribution to sharp resonances in the transmission at these 
cavity parameters. 

Inside the lead, we must distinguish between propagating and evanescent modes. The 
longitudinal component of the eigenstates in the leads, for propagating modes, can be written 




a 



•n 
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ik-n x 



(2.4) 
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where the wavevector, k n is given by k n = \J^r- — ( ! f : ) 2 - If there are N propagating modes 
then n — 1, N. Here we use a unit current normalization. A particle propagating in the 
n channel has energy in the interval, n 2 £i<£<(n + l) 2 £i, where E x = |^ = 0.0738eV. 
All remaining modes, n = N + 1, oo, are evanescent and can be represented in terms of 
the states, 

T nXn (x) = -^e—W, (2.5) 



where K n = \J (^) 2 — ^jjjf-- in Ref. [19], we showed that the evanescent modes for this system 
do not make a significant contribution to the scattering properties that we are interested in. 
Therefore we will neglect the evanescent modes in subsequent sections. 

As shown in Ref. [19], The eigenvalue equation, H\E) = E\E), leads to the relation, 



r„x„(o) = = E &..»' f % ) r »'' ( 2 - 7 ) 



Also, continuity of energy eigenstates at x = gives 

M 

j=l n'=l \ ax / x =o 

where 

**" - tag (B-A,) (2 ^ 8) 

is the (n, n') t/l matrix element of the reaction matrix [17]. The quantity, 0j,„(O) is a measure 
of the overlap between the j th cavity state, and the n th channel in the lead, evaluated at the 
interface, 

<M0) = yfjf dy 0,(0,1/) sin(^). (2.9) 

Let us now form an A^xl column matrix, b (a) whose matrix elements consist of the 
N probability amplitudes {b n } ({a n }) of the outgoing (incoming) propagating modes. The 
waveguide scattering matrix (S-matrix), S, is a NxN matrix which connects the incoming 
propagating modes to the outgoing propagating modes, b = S-a. 
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For the case when there are N propagating modes in the lead we can obtain an NxN 
S-matrix which may be written 

- _ (In - iK) 

S ~ (lN + iKY (2 ' 10) 
where Ijv is NxN unit matrix, and the NxN matrix, K, has matrix elements, K njJl / = 
\/knRn,n'\fkn' an d can be written 

K = wt- 1 - -w. (2.11) 

hii-M — tl in 

In Eq. (2.11), 1m is the MxM unit matrix, Hi n is an MxM diagonal matrix formed with 
the eigenvalues, Xj (j — 1, ...M) in the cavity, and w isan MxJV coupling matrix, 

/ Wi tl . . . W, N \ 



w = 



(2.12) 



••• W M ,N J 

where w^ n > = (f>j, n '(^)VKJ- With some algebra, the S-matrix can also be written in the form 

S = - (l N - 2iw^—= -r r w\ . (2.13) 

V E\ M - H m + iw-w^ J v ; 

In Ref. [19], we showed that, if evanescent modes are included, an additional term appears 
in the denominator in Eq. (2.13). 

The reaction matrix approach to waveguide scattering provides a very efficient means 
of computing the statistical properties of the scattering process because the Schrodinger 
equation only needs to be solved once to obtain the basis states and eigenvalues in the 
cavity. Using these values, the S-matrix can then be obtained at all other particle energies, 
E. Typically for the ripple cavity in Figure (1), we can obtain the scattering matrix at 10 5 
different values of incident energy in a reasonable amount of time on a Cray machine. 

One of the goals of this paper is to compare the statistical properties of the deterministic 
scattering process in the ripple cavity to statistical properties of a hypothetical scattering 
process in which H in is replaced by a diagonal matrix, H goe , composed of the M eigenvalues 
of an MxM Gaussian Orthogonal Ensemble (GOE) Hamiltonian, H' goe , and the iV columns 
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of the coupling matrix, w, are replaced by iV of the M normalized eigenvectors of H' to 
yield a coupling matrix, w goe [32,4]. In this random matrix theory approach, the strength of 
the coupling between the cavity and the lead is given by the phenomenological parameter, g. 
The parameter, g, does not appear in the deterministic scattering process. For deterministic 
scattering the strength of the coupling is entirely determined by the geometry and the 
potentials at the interface. The scattering matrix, obtain from RMT, can then be written 



>goe 



-(l N - 2igwl~ 1 — -w goe ). (2.14) 

V El M - H goe + igW goe -Wgoe / 



It was shown in Ref. ( [28]), using supersymmetry techniques, that for the case when the 
distribution of energy eigenvalues of H goe is centered at E — and M^oo, the average 
S-matrix can be written 



{Sgoe} = S1 N With S = - ; J; - /o ; _, , /r7nl i (2.15) 



1 - g[iE /2 + irv(E)} 
1 + g[iE/2 + Trv(E)] ' 

where v(E) = — (E/2) 2 is the average density of energy eigenstates. It is useful to 

introduce the quantity 

1 - s E . . 

H = fi r + iHi= ^ = gun + ig — , (2.16) 

where [i r = gvn and /ij = g^, respectively, are the real and imaginary parts of fi. The case 
when g = 1, corresponds to ideal coupling. In the neighborhood of £^0, the eigenvalues 
of H goe have a constant density, j-, and the average S-matrix (S goe ) = 0. When g^l, the 
average S-matrix cannot be zero. 



III. EIGENPHASES OF THE SCATTERING MATRIX 

We have analyzed some of the statistical properties of the eigenphases of the S-matrix 
for the case of deterministic scattering from the ripple cavity for the cases when the internal 
dynamics in the cavity is completely chaotic and when it is near-integrable. In this section, 
we consider the energy interval 256Ei<E<289E 1 when 16 channels are present in the lead. 
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The S-matrix is a 16x16 matrix, and for each value of incident energy it has 16 eigenvalues, 
e l5a (a = 1, 16) and 16 eigenfunctions, \5 a ) (S\5 a ) = e i5< * |£q,)). The S-matrix is unitary so 
the eigenfunctions, \5 a ), form a complete orthonormal set. We can use the orthonormality 
of the eigenfunctions to follow each eigenphase, 5 a , continuously as a function of energy [29]. 
The eigenfunctions, for two S-matrices evaluated at nearby energies, will be approximately 
orthogonal if they do not belong to the same eigenphase. Thus we can plot each of the 16 
different eigenphases as a function of the incident energy. These are shown in Figure (3) 
where the eigenphases, which are defined mod 2%, are "unwrapped" and allowed to evolve 
continuously as a function of energy. In Fig. (3. a), we show the case with ripple amplitude, 
a = 25A, where the classical cavity dynamics is chaotic, and in Fig. (3.b) we show the 
case a = lA where the classical cavity dynamics is mixed (see Fig. (2)). The case of mixed 
dynamics shows many more abrupt changes of phase as a function of energy than the chaotic 
case. This is due to the fact that the mixed dynamics has many long lived resonances not 
found in the chaotic case. This was also seen in Ref. [29,26]. We shall return to this feature 
in Section (IV). 

Below we first discuss the effect of direct processes on the distribution of eigenphases, 
and then we determine the distribution of nearest neighbor spacings of these eigenphases. 



A. Distribution of Eigenphases 

When a scattering process has a non-zero average S-matrix, (S), it indicates that direct 
processes may play a significant role in the scattering process. Direct processes are generally 
scattering events which do not interact significantly with the reaction region (cavity) [34]. 
When direct processes are present, the distribution of S-matrix elements that minimizes 
information about the scattering process is the Poisson kernel. For the case of an N channel 
process whose dynamics is time reversal invariant, the Poisson kernel has the form, 

P f ._l MlOT)g (3D 
^ N[b) ~Q \Det(l-(S)*S)\^) ' [6 ' 
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where fiis a normalization factor that ensures that the Poisson kernel satisfies the normal- 
ization condition, JcISPn(S) = 1. 

The S-matrix can be diagonalized by a unitary matrix, U, and, as mentioned earlier, the 
eigenvalues of the S-matrix are denoted, e l5a , a — 1, N. In terms of the eigenphases, 5 a , 
the normalization condition for the Poisson kernel, Eq. (3.1), can be written 



where Qjj is a normalization constant. Note that Pn(5i, 5n) is the joint probability 
density to find the angles, 5 a , in the intervals S a — >5 a + d5 a , (a = 1, N). 
In Ref. [37], it is shown that if the following change of angles is introduced 



which is just the distribution for the Circular Orthogonal Ensemble (COE) [5], [6]. Thus, 
even for scattering processes which include direct processes, it is possible in principle to 
transform away the direct processes and compare the eigenphase distribution with that of 
COE (note that similar ideas first appeared in literature in Ref. [15]). It is important to 
note that the transformation that removes direct processes is not the same as the unfolding 
process that occurs on the energy spectrum of bounded systems to give a constant average 
density. 

In our subsequent analysis, the case of scattering with only a single channel will be 
useful for analyzing data. For single channel scattering (TV = 1), the S-matrix reduces to 
the complex function, S = e lS , and the Poisson kernel, reduces to [14] 




(3.2) 




(3.3) 



and if one assumes ideal coupling, g — 1, then Eq. (3.2) reduces to 




(3.4) 
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Pis{6) 



1 [(1 - s*s] 
2^|(l-s*e i5 )| 2 ' 



(3.5) 



with normalization condition, J* n d5Pi(5) = 1. Under the transformation above, 
Pi(<J)->Pi(0) = 1/2tt, which is the COE prediction. 

Having obtained numerical values of the eigenvalues, e tSa , as a function of energy, we can 
compute an average value for each of the 16 eigenvalues, 



where 77 is the number of energy values used. The apparently continuous eigenphase curves 
actually consist of about 40, 000 discrete energy points. The approximate orthonormality of 
S-matrix eigenvectors for neighboring energies has been used to sort the eigenphases. Thus, 
the eigenphases and eigenvectors have an energy interval over which they are correlated, 
and we have used that fact in our sorting process. On the other hand, this correlation of 
the S-matrices at neighboring energies can prevent us from obtaining statistics that can be 
compared to RMT predictions. Comparison to RMT requires use of independent data points. 
Therefore, in order to study the statistical properties of the eigenphases, we must choose 
values of the eigenphases separated in energy a distance greater than the correlation length. 
For each eigenphase curve we select points which have an energy spacing, AE = 0.495£'i. We 
choose this spacing based on an analysis of the delay time correlation discussed in Section 
IV. (The delay time auto-correlation function is the second derivative of the eigenphase 
auto-correlation function.) 

We have computed a histogram of the number of eigenphases, N(S) = t)Pn(S), versus 
value of eigenphase, S a , in the 16 channel region, where 77 is the number of data points. 
In order to improve the statistics, we use data from four different ripple amplitudes, a = 
25A, 30A, 35A, 45A, all of which lie in the chaotic regime. All eigenphases lie in the 
energy interval, 256 < E/E 1 < 289 and have energy spacing, AE = 0.495£'i. Thus the 
histogram includes 67X16X4 = 4288 data points. We have found that the distribution of 
eigenphases, along a given eigenphase curve, is well described by the Poisson kernel for the 



1 v 



(3.6) 
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single channel case. We proceed as follows. We compute the average eigenvalue, s a = (e t5a ), 
for each eigenphase curve. We can form a histogram using the 67 data points from a single 
eigenphase curve. We do this for each of the 64 eigenphase curves, using the same number 
of bins and bin width for each curve. We then add these 64 histograms together to form 
a single histogram, which is shown in Fig. (4. a). The solid line in Fig. (4. a), is the single 
channel Poisson kernel, -Pi (5), but with (s) = ^Ea=i s «- We can use Eq. (3.3) to transform 
away the effects of the direct interactions and find the distribution of the transformed angles, 
9 a . If we use Eq. (2.16) and obtain \x from the numerically calculated values of s a = (e l5a ), 
then the transformed angles, 9 a , are given by 



with fi r and fa computed numerically from s a . In Fig. (4.b), we show how the histogram 
in Fig. (4. a) changes if we transform each eigenphase, 5 a , using Eq. (3.7). In this case 
distribution is approximately constant and equal to the area, Thus, having transformed 
away the contribution from the direct processes, we obtain the COE eigenphase distribution 
for this chaotic scattering process, with fairly high confidence level. 

For a = 1A the plots of the eigenphase distributions look very similar to the chaotic case 
shown in Fig. (4), and it appears that the eigenphase distribution is not as sensitive an 
indicator for underlying chaos as is the nearest neighbor spacing distribution, at least with 
this type of analysis. 



In this section we consider the nearest neighbor spacings between eigenphases, 5 a , of 
the scattering matrix for the 16 channel case in the energy interval 256 < E/Ei < 289. 
For any given value of the energy, the S-matrix only has 16 eigenphases. However, we can 
form a histogram of nearest neighbor eigenphase spacings if we obtain eigenphase spacings 
for a sequence of different energies in the range, 256 < E/E 1 < 289. In Fig. (5. a) we 




(3.7) 



B. Nearest Neighbor Eigenphase Spacing 
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show histograms of 1005 (15X67) nearest neighbor spacings for eigenphases computed at 
energy increments, AE = 0.495i?i and obtained by averaging over histograms for each ripple 
amplitude, a = 25A, 30A, 35A, 45A We have fit the histogram to the Brody distribution, 
[32] 

/ a \P r / a r 1 /2 + 3\-\ 1+ P 

w=i, y ap i- < y i with {= Miri)] ■ (3 ' 8> 

(a) is the average spacing between nearest neighbor eigenphases, and T(x) is the Gamma 
function. In Fig. (5. a), the solid line is a fit to the Brody distribution for (3 = 0.635. In 
Fig. (5.b) same calculation is performed after direct processes were transformed away. In this 
case P = 0.865. Note that the GOE prediction ((3 = 0.95) for the closed system eigenvalue 
spacings is fairly close to our value of (3, after the effects of direct processes are transformed 
away. In Ref . [30] , the nearest neighbor energy eigenvalue spacings for a closed ripple billiard 
were fit to the Brody distribution with f3 = 0.806. In that case, the deviation from GOE 
predictions was found to be due to bouncing ball orbits. Our result also contains bouncing 
ball contributions. It is useful to note that in our scattering system there is no long range 
energy correlation for 16 channel region (we have explicitly removed energy correlations by 
taking data points at large energy increments), in contrast to the case reported in reference 
[31] 

We also obtained a nearest neighbor spacing histogram for the case with mixed phase 
space in the 16 channel energy interval. In Fig. (6. a) we show the histogram of nearest 
neighbor eigenphase spacings for a = 0.5A, lA, 2A, 3A, 4A, 5A in the energy interval 
256 < E/Ei < 289 before direct process are transformed away. We use an energy spacing 
AE = 0.495.&L and obtain 1005 data points for each amplitude. We then average over the 
histograms for the six ripple amplitudes. The solid line is a fit to the Brody distribution for 
(3 = 0.2. In Fig. (6.b), the same histogram is shown after direct processes are transformed 
away. The Brody parameter in this case, (3 = 0.116, which is closer to a Poisson distribution 
(the Brody distribution becomes a Poisson distribution for (3 = 0). 

The distribution of nearest neighbor eigenphase spacings has been computed for a energy 
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independent scattering matrix in [39] for a very different physical system. They also report 
close agreement with COE predictions for the chaotic region, although they do not have to 
deal with direct processes. It is clear that direct processes can play an important role in 
causing deviations from random matrix theory predictions for scattering processes. 

IV. PARTIAL DELAY TIMES 

In this section we compare the partial delay time distribution, computed for the de- 
terministic scattering process, to values obtained from random matrix theory. The partial 
delay times are given by the energy derivative of eigenphases, r a = h%t [35]. The average 
partial delay time density for a scattering process governed by the scattering matrix, S goe , 
has been computed by [36,37], using supersymmetry techniques, and is given by, 

nM h//\nV/^ t)\ (l/2) jV/2 exp(-l/2(r/(r))) 

P(r) = (l/N) 2_j (<Hr - r a )) goe = r(iV /2) ( r /( r ))JV/2+ 2 (41) 

where (r) = l/N. 

Before showing the distribution obtained for the partial delay times, it is useful to 
discuss energy correlations contained in the partial delay time curves. In Fig. 7, we 
show the auto-correlation function for the partial delay times obtained in the 16 chan- 
nel energy regime 256 < E/E x < 289 and averaged over 6 different ripple sizes, a = 
22A, 23A, 24A, 25A, 26A, 27 A. For each partial delay time curve we obtain a auto- 
correlation function, and then we average over all 96 curves. We also show the GOE predic- 
tion for the 16 channel case as well as the partial delay time auto-correlation function for 
the near integrable regime (a = I A). The GOE prediction is obtained after performing the 
triple integration given in reference [33]. The energy scale is adjusted to correspond to the 
relevant scale for our data. We also note the partial delay time auto-correlation function is 
the second derivative of the eigenphase auto-correlation function. Therefore the eigenphase 
auto-correlation function decays more slowly than the partial time delay auto-correlation 
function. 
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In Fig. (8. a), we show a histogram of the scaled partial delay times, r/(r). We again 
consider the energy regime with 16 channels and vary the energy in the interval, 256Ei < 
E < 289-Ei. To obtain enough values to build good statistics, we use 4 different ripple 
sizes, a = 25A, 30A, 35A, A5A. For these ripple amplitudes, the ripple cavity dynamics is 
chaotic. We used 100 energy points per specific ripple size, and therefore an energy increment 
of AE = 0.33i?i. The average, delay time, (r), is obtained numerically for each partial delay 
time curve. Then the histograms for 64 scaled partial delay times are combined into one 
histogram by simply adding values in the corresponding bins. The solid line in Fig. (8. a), 
is a plot of the RMT prediction, N(t) = rf p(j), where r( is the area under the curve 
and (t) = l/N. The agreement is not good because our data contains the effect of direct 
scattering processes. In Fig. (8.b) we show the partial delay time density obtained from the 
eigenphases, 9 a , which no longer contain the effect of direct scattering processes. The solid 
line is a plot of N(t) with (t) = l/N [36]. The agreement is very good. Finally in Fig. (8.c) 
we plot histogram of 4000 partial delay times obtained from a 16 x 16 S goe by using different 
realizations. Again, the agreement between the data and Eq. (4.1) is very good. Thus, after 
the removal of the effects of direct scattering processes our deterministic scattering from 
the chaotic ripple cavity behaves very much like the RMT prediction. (It is useful to note 
that in Ref. [19], we compared the Wigner- Smith delay time distribution with numerically 
computed predictions of RMT. The Wigner-Smith delay time is defined, t ws = jjJ2a=i T a-) 

In Fig. 9, we show the delay time distributions for the near integrable 1A. We 

have used energy increments, AE = 0.1-Ei justified from Fig. 7. The delay time distribution 
for the near integrable case deviates significantly from the random matrix result and the 
results for chaotic cavity shown in Fig. (8). 

Let us now return to the eigenphase curve in Fig. (3.b). We see that the curves for the 
near integrable case have a sequence of fairly abrupt large changes of phase. These are due 
to resonance structures that cause larger than average delays of the particle in the cavity. In 
Fig. (10), we plot the Wigner-Smith delay time (which is an average over all partial delay 
times) in the energy interval, 257.4<£'/£' 1 <259.5. This energy interval contains two of the 
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large phase changes in the eigenphase curves in Fig. (3.b). We see that each large phase 
change gives rise to a large peak in the delay time. The crosses in Fig. (10) give the energies 
of the cavity basis states, Xj, in that energy interval. There appears to be one cavity state 
which lies at each resonance energy. In Fig. (10), we have also plotted the configuration 
space distribution of four of the cavity eigenstates, two at resonance and two off resonance. 
In Figs. (2.b)-(2.e), we have shown Husimi plots of the quantum Poincare surface of section 
for each of these four states. The two states giving rise to the delay time resonance peaks lie 
in the dominant nonlinear resonance structures in the classical phase space. The quantum 
particle appears to tunnel into these dynamical resonance structures, and is delayed there 
for a considerable length of time. 

V. CONCLUSION 

We have analyzed the statistical properties of a scattering process in a waveguide with 
a cavity which allows a range of dynamics, including integrable, mixed, or chaotic. In 
this waveguide, direct processes also play an important role. The "ripple" cavity that we 
use has the special feature that it allows us to form a Hamiltonian matrix to describe the 
dynamics interior to the cavity. This, in turn, allows us to use the reaction matrix approach 
to scattering for our deterministic scattering process. The reaction matrix approach is one 
of the most efficient methods for obtaining the large amounts of data necessary to obtain 
good statistics. Until now, mesh based models (like the boundary element method, finite 
element method, or recursive Green's function method) were the main numerical methods to 
deal with scattering problems. However, these methods use an energy dependent boundary 
condition which makes it a formidable task to obtain solutions for very large numbers of 
energy points. The reaction matrix approach allows us to circumvent this problem. It is 
also useful to note that the reaction matrix approach has been used extensively to study 
properties of the complex poles of the S-matrix. This is discussed in some detail in [39-42]. 

We have obtained a number of results. We find that, in the near integrable regime, 



16 



nonlinear resonances in the classical phase space give rise to large eigenphase excursions and 
long delay times for quantum particles that can tunnel into these dynamical structures. 

We have focused much of our discussion on the energy regime in which sixteen chan- 
nels are open in the lead. We have been able to follow each eigenphase of the S-matrix 
continuously as a function of energy. We have examined the statistical properties of the 
scattering process by gathering data about each eigenphase at discrete energy intervals in 
the sixteen channel regime. This assumes a kind of "stationarity" as a function of energy, 
of the underlying scattering process. We have chosen the energy intervals so that our data 
points are statistically independent. 

We have shown, for the scattering system considered here, that the affect of direct 
processes on the eigenphase curves can be transformed away. We find that, for the case 
where the cavity dynamics is classically chaotic, a partial time delay density histogram, 
formed from all sixteen transformed eigenphase curves, agrees to 96% confidence level with 
a Brody distribution with Brody parameter, (3 = 0.87. Similar deviations from the GOE 
prediction of (5 — 0.95 have been seen in the nearest neighbor energy eigenvalue spacing 
distributions of closed ripple billiards [30] and in that case are caused by bouncing ball 
orbits. We expect the same mechanism is having an affect here. 

VI. ACKNOWLEDGEMENTS 

The authors wish to thank the Welch Foundation, Grant No.F-1051; NSF Grant INT- 
9602971; and DOE contract No.DE-FG03-94ER14405 for partial support of this work. We 
also thank the University of Texas at Austin High Performance Computing Center for use 
of their computer facilities. 



17 



REFERENCES 

[1] E.P. Wigner, Can. Math. Congr. Proc. (University of Toranto Press, Toranto, 1957). 

[2] J. Garg, J. Rainwater, J.S. Petersen, and W.W. Havens, Jr., Phys. Rev. 134 B985 
(1964). 

[3] T.A. Brody, J. Flores, J.B. French, PA. Mello, A. Pandey, and S.S.M. Wong, Rev. Mod. 
Phys. 53 385 (1981). 

[4] M.L. Mehta, Random Matrices and the Statistical Theory of Energy Levels 2nd Ed. 
(Academic Press, New York, 1991). 

[5] F.J. Dyson, J. Math Phys. 3 140 (1962); 3 157 (1962); 3 166 (1962); 3 344 (1962). 

[6] L.K. Hua Harmonic Analysis of Functions of Several Complex Variables in the Classical 
Domains (American Mathematical Society, Providence, 1963). 

[7] S.W. MacDonald and A.N. Kaufman, Phys. Rev. Lett. 42 1189 (1979). 

[8] O. Bohigas, M.J. Giannoni, C. Schmit, Phys. Rev. Lett. 52 1. (1984) 

[9] G. Casati and F. Valz-Gris, and I. Guarnieri, Lett. Nuovo Cimento 28 279 (1980). 

[10] M.V. Berry and M. Tabor, Proc. Roy. Soc. Lond. A356 375 (1977). 

[11] L.E. Reichl, The Transition to Chaos in Classical Conservative Systems: Quantum 
Manifestations (Springer- Verlag, New York, 1992). 

[12] T. Guhr, A. Muller-Groeling, and H.A. Weidenmuller, Phys. Rept. 299 189 (1998). 

[13] C. Mahaux and H.A. Weidenmuller Shell Model Approach to Nuclear Reactions (North- 
Holland, Amsterdam, 1969) 

[14] G. Lopez, P. A. Mello and T. H. Seligman, Z. Phys. A 302 351 (1981) 

[15] P. A. Mello, P. Pereyra and T.H. Seligman, Annals of Physics 161 254 (1985), W. A. 
Friedman and P. A. Mello, Annals of Physics 161 276 (1985) 

18 



[16] P.W. Brouwer, Phys. Rev. B 51 16878 (1995) 

[17] E.P. Wigner and L.E. Eisenbud, Phys. Rev. 72 29 (1949). 

[18] G.A. Luna-Acosta, Kyungsan Na, L.E. Reichl, and A. Erokhin, Phys. Rev. E (1996) 53 
3271 (1996). 

[19] Gursoy B. Akguc and L.E. Reichl, Phys. Rev. E. 64 056221 (2001). 

[20] CM. Marcus, A.J. Rimberg, R.M. Westervelt, P.F. Hopkins, and A.C. Gossard, Phys. 
Rev. Lett. 69 506 (1992). 

[21] A.M. Chang, H.U. Baranger, L.N. Pfeiffer, and K.W. West, Phys. Rev. Lett. 73 2111 
(1994). 

[22] J. P. Bird, J. Phys.: Condens. Matter 11 R413 (1999). 

[23] H-D. Graf, H.L. Harney, H. Lengeler, C.H. Lewenkopf, C. Rangacharyulu, A. Richter, 
P. Schardt, and H.A. Weidenmuller, Phys. Rev. Lett. 69 1296 (1992). 

[24] J. Stein, and H.-J. StSckmann, Phys. Rev. Lett. 68 2867 (1992). 

[25] H.-J. Stockmann, Quantum Chaos - An Introduction (Cambridge University Press, Cam- 
bridge, UK, (1999)). 

[26] G.A. Luna-Acosta, J. A. Mendez-Bermudez, P. Seba, and K. N. Pichugin, Phys. Rev. 
E (1996) 65 046605 (2002). 

[27] C. Jung, T. H. Seligman, J. Phys. A. 28 1507 (1995) 

[28] J.J.M. Verbaarschot, H.A. Weidenmuller, and M.R. Zirnbauer, Phys. Rep. 129 367 
(1985). 

[29] Gursoy B. Akguc and L.E. Reichl, J. Stat. Phys. 98 813 (2000). 

[30] Wenjun Li, L.E. Reichl, and Biao Wu, Phys. Rev. E 65 056220 (2002) 



19 



[31] B. Dietz, M. Lombardi, and T. H. Seligman J. Phys, A 29 L95 (1996) 
[32] T.A. Brody, Lett. Nuov. Cim. 7 482 (1973) 

[33] N. Lehman, D.V. Savin, V.V. Sokolov, H.-J. Sommers, Physica D 86 572 (1995) 

[34] P. A. Mello and H. U. Baranger Waves Random Media 9 105 (1999) 

[35] E.P Wigner, Phys. Rev. 98 145 (1955); F.T. Smith, Phys. Rev. 118 349 (1960); D. 
Bohm, Quantum Theory (Prentice-Hall, Inc. Englewood Cliffs, N.J., 1950), p. 260. 

[36] Y.V. Fyodorov and H-J Sommers, J. Math. Phys. 38 1918 (1997). 

[37] Dmitry V. Savin, Yan V. Fyodorov and H-J Sommers, Phys. Rev. E. 63 035202(R) 
(2001). 

[38] F. Haake, F. Izrailev, N. Lehman, D. Saher, and H.-J. Sommers. Z. Phys. B 88 359 
(1992) 

[39] C. Jung, M. Muller, and I. Rotter, Phys. Rev. E 60 114 (1999). 

[40] I. Rotter, E. Persson, K. Pichugin, and P. Seba, Phys. Rev. E 62 450 (2000). 

[41] I. Rotter, Phys. Rev. E 64 36213 (2001). 

[42] H.-J. Stockmann, E. Persson, Y.-H. Kim, M. Barth, U. Kuhl, and I. Rotter, Phys. Rev. 
E 65 66211 (2002). 



20 



FIGURES 

FIG. 1. The geometry of the two dimensional ballistic waveguide used in our calculations; a is 
the half-width of the ripple, d is the width of the lead and the average width of the cavity. The 
ripple cavity extends from x = to x = L 

FIG. 2. Surfaces of section for L = 51lA, d = lOlA, and a = lA. (a) A Poincare surface of 
section showing p x /p = cos(a) versus x, each time the particle hits the bottom wall, (b) Husimi 
plot of quantum surface of section (QSS) for cavity eigenstate with Xj = 257.1934i?i. (c) QSS for 
cavity eigenstate with Xj = 257.9339£i. (d) QSS for cavity eigenstate with Xj = 258.6655Ei. (e) 
QSS for cavity eigenstate with Xj = 258.9072£i. 

FIG. 3. Eigenphases, 5 a versus E/E\ for the energy interval, 256<£'/£i<272: (a) a = 25 A, 
(b) a = lA. 

FIG. 4. (a) Histogram of number of eigenphases, N(5), versus S for the 16 channel energy 
interval 256<E/E 1 <289, and for four different ripple sizes, a = 25A, 30 A, 35A, 45A. The solid 
line is a plot of the single channel Poisson kernel with (S) = X^=i(e Jl5< *) and normalized to the 
number of eigenphases. (b) Histogram of transformed eigenphases, 9 a , for all 64 eigenphase curves. 
For all cases, d = 101 A and L = 511 A. A \ 2 test result is also shown for both plots with 17 bins 
taken into account 

FIG. 5. Histogram of number, N(a), of nearest neighbor scaled eigenphase spacings, a, for the 
chaotic regime with d = lOlA, L = 51lA. The average spacing, (a) is obtained for each eigenphase 
curve. The histograms contain a total of 15x67 = 1005 data points averaged over four different 
ripple sizes, a = 25A, 30A, 35A, 45A. (We obtain a histogram for each of the four values of 
the ripple amplitude. We then add them and divide by four.) (a) Before direct processes are 
transformed away. The thin solid line is the Brody distribution with (3 = 0.635. (b) After direct 
processes transformed away. The thin solid line is the Brody distribution with [3 = 0.865. A x 2 test 
result is also shown for both plots with 13 bins in it taken into account 
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FIG. 6. Histogram of number, N(cr), of nearest neighbor eigenphase spacings, a, for the near 
integrable regime with d = lOlA, L = 511A The histograms contain a total of 1005 data points 
averaged over 6 different ripple sizes, a = 0.5A, lA, 2A, 3A, 4A, 5A. (a) Before direct processes 
are transformed away. The thin solid line is the Brody distribution with (3 = 0.2. (b) After direct 
processes are transformed away. The thin solid line is the Brody distribution with f3 = 0.116. A 
X 2 test result is also shown for both plots with 8 bins taken into account 

FIG. 7. The auto-correlation function of time delays in the 16 channel energy interval, 
256 < E/E\ < 289. The thin line is obtained numerically for the chaotic regime, using six dif- 
ferent ripple amplitudes, a = 22A, 23A, 24A, 25A, 26A, 27 A, with direct processes transformed 
out of the data. An auto-correlation function is obtained for each partial time delay curve and 
the average of those 96 auto-correlation functions is shown. The dotted-dashed line is the GOE 
result for perfect coupling with 16 modes with the average density of states chosen equal to 1.25. 
The thick line represents the numerically obtained auto-correlation function in the region of mixed 
phase space for a = 0.5A, 1A, 2A, 3A, 5A 

FIG. 8. Histogram of number of scaled partial delay times, N(t), versus r/(r), for the 16 
channel energy interval, 256<E/E\<289, with d = WlA and L = 511A Data for ripple amplitudes 
a = 25A, 30A, 35A, 45A is included in the histograms. Data points are taken at energy intervals, 
AE = 0.33£i. (a) Histogram of scaled partial delay times taken from eigenphase curves for 5 a . 
A scaling factor, (r), is obtained for each eigenphase curve, (b) Histogram of scaled partial delay 
times taken from eigenphase curves for the transformed eigenphases, 9 a . (c) Histogram of partial 
delay times obtained from the 16x16 S-matrix, S goe (includes 4000 data points). A \ 2 test result 
is also shown for the plots (b) and (c) with 13 bins taken into account. 
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FIG. 9. Histogram of partial delay times for the 16 channel energy interval, 256<-E/£i<289, 
with d = lOlA and L = 51lA. Data for ripple amplitudes, a = 0.5A, lA, 2A, 3A, 5A, is used to 
construct the histograms. Data points are taken at energy spacings, AE = 0.1E±. A total of 
400X16 points is used, (a) Histogram of scaled partial delay time curves taken from eigenphase 
curves for S a . (b) Histogram of scaled partial delay times taken from curves for transformed 
eigenphases, 9 a . 

FIG. 10. Plot of Wigner-Smith delay time in the energy the interval 257 A<E / E±<259. 5 for 
d = lOlA, L = 511A, a = lA. Crosses show values of cavity basis state energies in this interval. 
Inserts show the spatial distribution of four cavity basis states, two at resonance and two off 
resonance. 
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